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Abstract 

This paper deals with quantum fluctuations near the classical instanton configuration. Feynman 
diagrams in the instanton background are used for the calculation of the tunneling amplitude 
(the instanton density) in the three-loop order for quartic double-well potential. The result for the 
three-loop contribution coincides in six significant figures with one given long ago by J. Zinn-Justin. 
Unlike the two-loop contribution where all involved Feynman integrals are rational numbers, in 
the three-loop case Feynman diagrams can contain irrational contributions. 
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Introduction 


0 , 


Euclidean classical solutions of the field equa- 

0-P 


There is no question that instantons 

tions, represent one of the most beautiful phenomena in theoretical physics bt-lai- In¬ 
stantons in non-Abelian gauge theories of the QCD type are important component of the 
non-perturbative vacuum structure, in particular they break chiral symmetries and thus sig¬ 


nificantly contribute to the nucleon (and our) mass |4j. Instantons in supersymmetric gauge 
theories lead to derivation of the exact beta function B, and in the “Seiberg-Witten” J\f=2 
case to derivation of the super potential by the exact evaluation of the instanton contribu¬ 
tions to all orders j(3]. The instanton method now has applications in stochastic settings 
beyond quantum mechanics or field theories, and even physics — in chemistry and biology - 


see e.g. discussion of its usage in the problem of protein folding in 


Since the work by A. Polyakov a the problem of a double well potential (DWP) has been 
considered as the simplest quantum mechanical setting illustrating the role of instantons in 
more complicated quantum field theories. In the case of the DWP one can perform certain 
technical tasks - like we do below - which so far are out of reach in more complicated/realistic 
settings. 

Tunneling in quantum mechanical context has been studied extensively using WKB and 
other semiclassical means. The aim of this paper is not to increase accuracy on these 
quantum-mechanical results, but rather to develop tools - Feynman diagrams on top of an 
instanton - which can be used in the context of many dimensions and especially in Quantum 
Field Theories (QFTs). Therefore we do not use anything stemming from the Schrodinger 
equation in this work, in particularly do not use series resulting from recurrence relations 
or resurgence relations (in general, conjectured) by several authors. 

Another reason to study DWP is existing deep connections between the quantum mechan¬ 
ical instantons - via Schrodinger equation - with wider mathematical issues, of approximate 
solutions to differential equations, defined in terms of certain generalized series. A particular 
form of an exact cuantization condition was conjectured by J. Zinn-Justin and collaborators 
(for a review see 


0 . 


12] and references therein), which links series around the instantons with 


the usual perturbative series in the perturbative vacuum. Unfortunately, no rigorous proof 
of such a connection exist, and it remains unknown if it can or cannot be generalized to 
the field theory cases we are mainly interested in. Recently, for the quartic double well and 
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Sine-Gordon potentials Dunne and Unsal (see [8] and also references therein) have presented 
more arguments for this connection, which they call resurgent relation. 

In { 9 ] the method and key elements (a non-trivial instanton background and new effective 
vertices) to calculate the two-loop correction to the tunneling amplitude for the DWP were 
established. In particular, the anharmonic oscillator was considered in order to show how to 
apply Feynman diagrams technique. In 10] the Green function in the instanton background 
was corrected, and it was attempted to obtain two- and three-loop corrections. Finally, 
Wohler and Shuryak [ll] corrected some errors made in 10] and reported the exact result 
for the two-loop correction. 

The goal of the present paper is to evaluate the three-loop correction to the tunneling 
amplitude and compare it with the results obtained in 12] by a completely different method, 
not available in the held theory settings. 


Three-loop correction to the instanton density 


Let us consider the quantum-mechanical problem of a particle of mass m — 1 in a double 
well potential 

V = A (x 2 — rj 2 ) 2 . ( 1 ) 

The well-known instanton solution X mst (t) = 77 tanh(|o;(t — t c )), with cu 2 = 8 A 77 s , describ¬ 
ing the barrier tunneling is the path which possesses the minimal action So = S[X inst (t)] = 

3 

Setting co = 1, and shifting coordinate to the minimum of the potential one gets 
the anharmonic oscillator potential in a form V av h = \x 2 — V2~\x 3 + A x 4 with one (small) 
dimensionless parameter A. J. Zinn-Justin et al 12] use the same potential with A = g/ 2. 
The classical action So of the instanton solution is therefore large and i§ used in the 
expansion.The ground state energy Eq within the zero-instanton sector (pure perturbation 
theory) is written in the form 


1 “ A 

rr — 1 X An 

0 9 / y on 1 


2 ^ SZ 

n =0 u 


(A) — 1 ) , 


( 2 ) 


Another series to be discussed is the splitting 5 E = E first excited state - E ground state related 
to the so called instanton density 19] in the one-instanton approximation as 


R 

SE = A 

n=0 D 0 


(Bo = 1 ) , 


( 3 ) 
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where A E = 2 A 4 e~ s ° is the well-known one-loop semiclassical result Coeficients A n 
in the series 02]) can be calculated using the ordinary perturbation theory (see 15]) while 
many coefficients B n in the expansion (|2]) were found by J. Zinn-Justin, 1981-2005 (see 


12 ] and references therein), obtained via the so called exact Bohr-Sommerfeld quantization 


condition. 

Alternatively, using the Feynman diagrams technique Wohler and Shuryak ll| calculated 
the two-loop correction Bi = —71/72 in agreement with the result by J. Zinn-Justin 12]. 
Higher order coefficients B n in ([3]) can also be computed in this way. Since we calculate the 
energy difference, all Feynman diagrams in the instanton background (with the instanton- 
based vertices and the Green’s function) need to be accompanied by subtraction of the 
same diagrams for the anharmonic oscillator, without the instanton (see [oj] for details). For 
r 1 it permits to evaluate the ratio 

(~v\e~ HT \v) 


i inst 


d -Hi 


| V) anh 

where the matrix elements (— r]\e~ HT \rj) inst , {y\e~ HT \i]) an h are calculated using the 
instanton-based and vacuum diagrams, respectively. 

The instanton-based Green’s function G(x , y) 


G(x,y) = G°(x,y) 


2-xy+-\x-y\(ll-3xy) + (x - yY 


+-( 1-* 2 )(1 -y 2 ) 


log G°(x,y)- 
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(4) 

is expressed in variables x = tanh(h-), y = tanh(y) , in which the familiar oscillator Green 
function of the harmonic oscillator is 


G°(x,y) = 


1 — \x — y\ — xy 


(5) 


1 + \x — y\ — xy 

In its derivation there were two steps. One was to fold a function which satisfies the Green 
function equation, used via two independent solutions and standard Wronskian method. The 
second step is related to a zero mode: one can add a term 0 o(G)</>o(O) with any coefficient 
and still satisfy the equation. The coefficient is then fixed from orthogonality to the zero 
mode, see 1 ( 1 . 

n 

The two-loop coefficient is given by the two-loop diagrams (see Fig. [Hand [111])- 


B\ — a + b\ + 62 
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where 


97 


a = 


bi = 


53 


bo = 


1680 ’ 1260 ’ 
reflecting the contribution of four Feynman diagrams. 


39 

560 ’ 


c = — 


49 
60 ’ 


( 6 ) 


The three-loop correction B 2 (j3J) we are interested in is given by the sum of eighteen 3-loop 
Feynman diagrams, which we group as follows (see Figs. [2]-|3]) 


(7) 


B 2 = Ui + fen + b\ 2 + b 2 \ + b 22 + 623 + ^24 

+ d-\-e + f + g + h +Ci+ C 2 + C 3 + C 4 + C 5 + C 6 + B 2 i oop , 
complementing by a contribution from two-loop Feynman diagrams 


„ 1, , , . 39589 

B 2100 P — - ( a + bi + b 2 ) + (a + 61 + b 2 ) c — 9 ^g 2 QQ > 

(see (jHJ)). Thus, all Feynman diagrams contributing to (J7]) are presented in Figs. dH The 
rules of constructing the integrals for each diagram should be clear from next two examples. 
The explicit expression for the Feynman integral b 23 in Fig. [2] is 
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»i 


623 = 8 L* 


dy / dz 


dw 


-1 






J{x,y,z,w) (xy zw G xx G xy Gy Z G yw G 2 zw — G 0 xx G 0 xy G 0 yz G 0 yw (G° zw ) 
while for C 4 in Fig. [3] it takes the form 

xy 




dy / dz- 


s~i s~i 

y2 \ ^xy ^yz ^ZZ 5 


( 8 ) 


(9) 


'-1 J-1 

here we introduced notations G xy = G(x,y), G® y = G°(x,y) and J = 
(1 _ 1 a ,2 ) ( jZyZ\ ( iz ~ z 2 ) 7 y ~2 , ■ Notice that the c's diagrams come from the Jacobian of the zero 
mode and have no analogs in the anharmonic oscillator (single well potential) problem. 

For calculation of the symmetry factors for a given Feynman’s diagram we use the Wick’s 


theorem and contractions, see e.g 


17]. It can be illustrated by the next two examples. For 


diagram b n the four propagators can be rearranged in 4! ways and the effect is duplicated by 
interchanging the two vertices t\, t 2 giving a symmetry factor of 2 x 4! =48. For the diagram 
C 4 the last propagator, which starts and ends at the same vertex forming loop, contributes 
with a factor of two which is also duplicated by rearranging the two inner propagators, 
finally, the symmetry factor of 2 x 2! = 4. 


5 









Results 


The obtained results are summarized in Table [U All diagrams are of the form of two- 
dimensional, three-dimensional and four-dimensional integrals. In particular, the diagrams 
b \i and d (see Fig. [2]) 


6ll = -r 

d = 


1 

48 

1 

16 


dx 


dy 


'-i 

r 1 


dx 


-i (1 -x 2 )(l - y 2 ) 
dy 


G 


xy 


{G\ 


o \ 4 

xy) 


( 10 ) 


( G G 2 G /hiu (G® i 2 g® \ 

l (j xx'^ X y^ jr yy ' jr xx\ {j *xy) Lz yy ) i 


f ^ J l — x 2 )[l — y 2 ) K Xy w \ xyj yy 

is given by two-dimensional integrals are the only ones which we are able to calculate ana¬ 
lytically 

1842223 I - h4l( 367 C(2) - 180 C(3) - 486 C(4)) = + b\\ rat 

J (11) 


6n = 
d = 


592704000 9800 

205441 525 


C(2) =d rat + d 


jirrat 


2469600 411600 

here ((n) denotes the Riemann zeta function of argument n (see E. Whittaker and G. Watson 
(1927)). They contain a rational and an irrational contribution such that 

^jirrat 


jirrat 

idrat 


-4.55 


d rat 


0.025 . 


It shows that the irrational contribution to bn is dominant with respect to the rational part 
while for diagram d the situation is opposite. Other diagrams, see Table I, were evaluated 
numerically with an absolute accuracy ~ 10 -7 . Surprisingly, almost all of them (15 out of 18 
diagrams) are of the same order 10 _1 as B 2 with few of them (diagrams ai, &12, &21) which 
are of order 10~ 2 . 

J. Zinn-Justin (see 12] and references therein) reports a value of 

6299 


r> Zinn—Justin 
^2 


10368 


-0.6075424 , 


( 12 ) 


while present calculation shows that 


££present 


-0.6075425 , 


(13) 


which is in agreement, up to the precision employed in the numerical integration. 

Similarly to the two-loop correction Bi the coefficient B 2 is negative. Note also that for 
Bi all diagrams are negative while for B 2 there are diagrams of both signs. For not-so- 
large barriers (So ~ 1), the two-loop and three-loop corrections are of the same order of 
magnitude. 
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The dominant contribution comes from the sum of the four-vertex diagrams 
£> 12 , & 21 , b 23 , e, h, ci, C 5 , C 6 while the three-vertex diagrams ai, b 22 , £> 24 , f, g, C 2 , C 3 , C 4 provides 
minor contribution, their sum represents less than 3% of the total correction i? 2 - Inter¬ 
esting that for both two and three loop cases the largest contribution comes from diagrams 
stemming from the Jacobian, c for Bi and C 5 , Cq for B 2 . Those diagrams are absent in the 
perturbative vacuum series, and thus do not have subtractions. 

We already noted that some individual three-loop diagrams contain irrational numbers. 
Since the J. Zinn-Justins result is a rational number, then there must be a cancelation of 
these irrational contributions in the sum (J7J). From (fill) we note that the term (b l {{ at + d irrat ) 
gives a contribution of order 10~ 2 to the mentioned sum (I7|) . and therefore the coincidence 
10 — ' between present result (fT31) and one of J. Zinn-Justin is an indication that such a 
cancelation occurs. Now, we evaluate the coefficients A\, A 2 in ([2]) using Feynman diagrams 
(see [3]). In order to do it let us consider the anharmonic oscillator potential V an h = 
\x 2 — y/2~\x 3 + Ax 4 and calculate the transition amplitude (x = 0\e~ HanhT \x = 0). All 
involved Feynman integrals can be evaluated analytically. In the limit r —> 00 the coefficients 
of order Sq 1 and Sq 2 in front of r gives us the value of A\ and A 2 , respectively. As it was 
mentioned above the c's diagrams do not exist for the anharmonic oscillator problem. The 
Feynman integrals in Fig. [Ogive us the value of A±, explicitly they are equal to 


1 



bi 


1 

24 ’ 



The diagrams in Fig. [2] determine A 2 and corresponding values are presented in Table Q] 
bn = — and d = — Straightforward evaluation gives 


Ai = 



1 

4 ’ 


which is in agreement with the results obtained in standard multiplicative perturbation 
theory (see 16]). No irrational numbers appear in the evaluation of Ai and A 2 . It is worth 
noting that (see Table |T]) some Feynman integrals being different give the same contribution, 


f = 9 = 


3 

32 


&22 — b 2 4 — 


1 

24 ' 


In the instanton background the corresponding values of these diagrams do not coincide but 
are very close. 
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Conclusions and Discussion 


In conclusion, we have calculated the tunneling amplitude (level splitting or the instanton 
density) up to three-loops using Feynman diagrams for quantum perturbations on top of the 
instanton. Our result for B 2 is found to be in good agreement with the resurgent relation 
between perturbative and instanton series suggested by J. Zinn- Justin (for modern reference 


see 


12 ]). 


Let us remind again, that this paper is methodical in nature, and its task was to develop 
tools to calculate tunneling phenomena in multidimensional QM or QFT context, in which 
any results stemming from the Schrodinger equation are not available. We use a quantum 
mechanical example as a test of the tools we use: but the tools themselves are expected to 
work in much wider context. 

One comment on the results is that the final three-loop answer has a rational value. How¬ 
ever, unlike the evaluation of the two-loop coefficient B\ where all Feynman diagrams turned 
out to be rational numbers, in our case of B 2 at least two diagrams contain irrational parts. 
What is the origin of these terms and how they cancel out among themselves are questions 
left unanswered above, since several diagrams had resisted our efforts to get the analytic 
answer, so that we used numerical multidimensional integration methods, in particular, a 


dynamical partitioning 


18]. Perhaps, this can still be improved. 


Another intriguing issue is the conjectured relation between the instanton and vacuum 
series: at the moment we do not understand its origin from the path integral settings. Some 
diagrams are similar, but expressions quite different and unrelated. New (tadpole) diagrams 
originate from the instanton zero mode Jacobian, and those have no analogues in the vacuum. 
Surprisingly, they provide the dominant contribution to two-, three-loop corrections B\ (one 
diagram out of four, see Fig.l) and B 2 (six diagrams out of 18, see Fig.3)): ~ 83% and ~ 97% 


20], respectively, see Table I. It implies that the sum of vacuum-like diagrams originating 


from the action is small. 

Finally, we note that to our knowledge this is the first three-loop calculation on a nontrivial 
background of an instanton. Similar calculations for gauge theories would be certainly 
possible and are of obvious interest. One technical issue to be solved is gauge Green function 
orthogonal to all (including gauge change) zero modes. 
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Note added in proof (July, 2015): 

During the time after the paper was submitted for publication we obtained a number of 
new results. We evaluated the contributions of c, C 5 -like diagrams, with maximal number 
of integrations, to the next order coefficients. The trend continues: those diagrams still 
contribute a significant fraction of the total answer, namely 83%, 127%, 60%, 20% of total 
two-, three-, four-, five-loop Bi, B 2 , B 3 , R 4 contributions, respectively. At the same time, 
surprisingly, the absolute values of all these diagrams are rather close. Advance in numerical 
multidimensional integrations lead to an increase in accuracy, the agreement in B 2 is now 
improved to six significant digits. 


Appendix (as of August 9, 2015) 


In this paper for the historical reasons we use the instanton-based Green’s function 


G(x,y) = G°(x, y) 


1 

2-xy+-\x-y\(ll-3xy) + (x - y) x 


+-(1-* 2 )(1 -y 2 ) 


log G°(x,y)-^- 


expressed in variables x = tanh(b-), y = tanh(y) , in which the oscillator Green function 


G°(x,y ) = 


1 — \x — y\ — xy 


1 + \x — y\ — xy ’ 

see ©. 

For the instanton field we have used the effective triple and quartic coupling constants 

v 3 = ~ tanh(«/2) V /2 U = \ So -1 , 

respectively ( S 0 is the action evaluated in the instanton solution, see p.3), while for the 
(subtracted) anharmouic oscillator we have 


R* = 


■£sr 


- 1/2 
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For tadpole c-diagrams the vertex is effectively represented by 


_ y/3 tanh(t/2) i/ 2 
tad ~ 4 cosh 1 2 3 4 5 (t/2) 0 ‘ 

The normalization of the Green function (3j) was chosen in such a way that in the r.h.s. 
of the defining equation contains the coefficient 4 in front of the delta function, |<5(x — y). 
Vertices V 3 , V 4 , V tad we have chosen accordingly. 

Certainly, the instanton-based Green’s function can be derived for the standard normaliza¬ 
tion with the coefficient 1 in front of the delta function S(x — y) in the defining equation. 
This instanton-based Green’s function corresponds to 


G(x,y) -A G°(x,y) 


2-xy + -\x — y\(ll — 3xy) + (x — y)“ 


+ -(l-s 2 )(l~y 2 ) 


l°g(2 G°(x, y)) ■ 
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where the oscillator Green function changes to 

1 — \x — y\ — xy 


Gr(x,y) -A 


2 (1 + |x — y\ — xy) 


= G„(x,y) . 


= G n (x, y) , 

(14) 

(15) 


Vortices have to be replaced as well 


V 3 ^2V2V 3 = V 3 n , V 4 ^4V 4 = V 4 n , V tad ^V2V tad = V t 


n 

tad 


(16) 


It can be easily checked that any Feynman integral is invariant under the transformations 
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Figure 1: Diagrams contributing to the two-loop correction B\ = a + b\ + 62 + c. They enter 
into the coefficient B 2 via the term h^Zoop- For the instanton field the effective triple and quartic 
coupling constants (vertices) are V 3 = — ^ tanh(i/2) S 0 1 ^ 2 and V 4 = ^ Sg 1 , respectively, while for 
the (subtracted) anharmonic oscillator we have V 3 = — ^ S 0 1 ^ 2 and V 4 = ^Sg 1 , all marked by 
(filled) bullets. The tadpole in diagram c, which comes from the zero-mode Jacobian rather than 
from the action, is effectively represented by the vertex (Jacobian source) V ta d = cosilh// 9 ) 
marked (unfilled) open bullet. The signs of contributions and symmetry factors are indicated. 
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Feynman 

diagram 

Instanton 

b 2 

Vacuum 

A 2 

ai 

-0.06495185 

5 

192 

bi2 

0.02568743 

1 

64 

b>2\ 

0.04964284 

11 

384 

b>22 

-0.13232566 

1 

24 

b>23 

0.28073249 

1 

8 

^24 

-0.12711935 

1 

24 

e 

0.39502676 

9 

64 

/ 

-0.35244758 

3 

32 

9 

-0.39640691 

3 

32 

h 

0.31424977 

3 

32 

Cl 

-0.3268200 

- 

C2 

0.63329511 

- 

C3 

0.12657122 

- 

C 4 

0.29747446 

— 

C5 

-0.77100484 

— 

C6 

-0.80821157 

- 

hD 

0.0963 

7 

384 

I 3 D 

-0.0158 

19 

64 

hD 

-0.8408 

155 

384 


Table I: Contribution of each diagram in Fig 12] - [3] for the three-loop corrections B 2 (left) and A 2 
(right) with symmetry factor included. We write B 2 = ( B 2100 P+I 2 D+I 3 D+I 4 D ) where I 2 D, I 3 D, hD 
denote the sum of two-dimensional, three-dimensional and four-dimensional integrals, respectively. 
Similarly, A 2 = I 2 D + I 3 D + hD- The term B 2 ioop = 39589/259200 ~ 0.152735 (see text). 
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1 1 1 
16 16 8 



Figure 2: Diagrams contributing to the coefficient E> 2 - Triple and quartic vertices V 3 , V 4 are marked 
by (filled) bullets. The signs of contributions and symmetry factors are indicated. 
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Figure 3: Tadpole diagrams contributing to the coefficient B 2 . They come from the Jacobian of 
the zero mode and have no analogs in the anharmonic oscillator problem. Tadpole vertex Vt a d 
(Jacobian source) is marked by (unfilled) open bullet. The signs of contributions and symmetry 
factors are indicated. 
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